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^ ■ Abstract 
in 

■ A coupled-channel analysis of the ^^'^'^'^^0(p,y) data has been performed to determine the neu- 

tron transition strengths of 2f states in Oxygen targets, using the microscopic optical potential 
and inelastic form factor calculated in the folding model. A complex density- and isospin depen- 
dent version of the CDM3Y6 interaction was constructed, based on the Brueckner-Hatree-Fock 
calculation of nuclear matter, for the folding model input. Given an accurate isovector density 
I dependence of the CDM3Y6 interaction, the isoscalar {6o) and isovector {6i) deformation lengths 



of 2f states in 18-20, 22q j-^g^yg been extracted from the folding model analysis of the {p,p') data. A 
specific A^-dependence of Sq and 6i has been established which can be linked to the neutron shell 
closure occurring at N approaching 16. The strongest isovector deformation was found for 2f state 
' in ^'^O, with Si about 2.5 times larger than 6q, which indicates a strong core polarization by the 

• valence neutrons in ^'^O. The ratios of the neutron/proton transition matrix elements (Mn/Mp) 

^ . determined for 2f states in 18,20q j^g^yg bggn compared to those deduced from the mirror symmetry, 

using the measured B{E2) values of 2'^ states in the proton rich "'^'^Ne and ^'^Mg nuclei, to discuss 
the isospin impurity in the 2^ excitation of the A = 18, T = 1 and A = 20, T = 2 isobars. 
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I. INTRODUCTION 



Inelastic proton scattering has been used over decades as a very efficient tool to yield 
the nuclear structure information. In difference from the electromagnetic probes, protons 
interact strongly with both protons and neutrons in the target nucleus, and the neutron 
and proton transition strengths of a nuclear excitation could be reliably deduced from the 
{p,p') measurement, in terms of the neutron M„ and Mp matrix elements [1]. The knowledge 
of Mn and Mp can shed light into the relative contributions by the valence nucleons and 
the core to the nuclear excitation, and hence, provides important information on the core 
polarization by the valence nucleons which could eventually lead to changes in the shell 
structure. This topic has recently become of significant interest in the experimental studies 
with radioactive beams where the inelastic proton scattering can be accurately measured, in 
the inverse kinematics, for the short-lived unstable nuclei (see, e.g., Refs. [2, 3, 4, 5, 6] for the 
ip,p') measurements with the unstable Oxygen isotopes). With large neutron (or proton) 
excess in the unstable nuclei, such proton scattering data provide also a vital information 
for studying the isospin effects in the proton-nucleus interaction. Although the isospin 
dependence of the nucleon optical potential (OP), known by now as Lane potential [7], has 
been studied since a long time, few attempts were made to study the isospin dependence of 
the transition potential or form factor (FF) for inelastic scattering. As neutron and proton 
contribute differently to the nuclear excitation, the inelastic scattering FF contains also an 
isospin dependence [8] which determines the degree of the isovector mixing in the inelastic 
scattering channel that induces the excitation [12] . 

In general, the isospin-dependent part of the nucleon-nucleus OP is proportional to the 
product of the projectile and target isospins, and the total OP can be written in terms of 
the isoscalar (IS) and isovector (IV) components [7] as 

U{R) = Uo{R) ± eUi{R), e = {N - Z)/A, (1) 

where the + sign pertains to incident neutron and - sign to incident proton. The strength 
of the Lane potential Ui is known from {p,p) and (n, n) elastic scattering and {p,n) reac- 
tions studies, to be around 30-40% of the Uq strength. In the first order of the collective 
model, inelastic nucleon-nucleus scattering cross section can be reasonably described, in the 
distorted-wave Born approximation (DWBA) or coupled channel (CC) formalism, with the 
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inelastic form factor F given by 'deforming' the optical potential (1) as 

The explicit knowledge of the deformation lengths and 5i would give us important struc- 
ture information about the IS and IV transition strengths of the nuclear excitation under 
study. There are only two types of experiment that might allow one to determine 5q and 5\ 
using prescription (2): 

i) Charge exchange (p, n) reaction leading to the excited isobar analog state. It was 
shown, however, that the calculated inelastic (p, n) cross sections were insensitive to 5\ due 
to the dominance of two-step process [9]. 

ii) Another way is to extract 5o(i) from the (p, p') and (n, n') data measured at about the 
same energy for the same excited state of the target [8, 9]. Since eUi/Uq is only about few 
percent, the uncertainty of this method can be quite large. Moreover, it remains techni- 
cally not feasible to perform simultaneously (p, p') and (n, n') measurements in the inverse 
kinematics for unstable nuclei. 

From a theoretical point of view, the form factor (2) has been shown to have inaccu- 
rate radial shape which tends to underestimate the transition strength, especially, for high- 
multipole excitations induced by inelastic heavy-ion scattering [10, 11]. As an alternative, a 
compact approach based on the folding model has been suggested in Ref. [12] for the deter- 
mination of the IS and IV transition strengths of the {AS — AT — 0) nuclear excitations 
induced by inelastic proton scattering. This compact folding approach was used with some 
success in the DWBA analyses of the ^^~'^^S{p,p') and ^^~^^0(p,p') data measured in the in- 
verse kinematics, to determine 6o and 6i for the 2^ states in the Sulfur and Oxygen isotopes 
under study [12, 13]. We recall that the basic inputs for such a folding -|- DWBA analysis 
are the effective NN interaction between the incident proton and those bound in the target, 
and the transition densities of the nuclear excitation. Consequently, for a carefully chosen 
model of the nuclear transition densities, the more accurate the choice of the effective A'^A'^ 
interaction the more reliable the deduced 6o and 6i deformation lengths. 

Our first folding model analysis of the ^^'^°0(p,p') data [12] has used a well-tested 
CDM3Y6 density dependent interaction [14] to construct the proton-nucleus OP and in- 
elastic FF. For simplicity, the density dependence of the isovector part of the CDM3Y6 
interaction has been assumed in Ref. [12] , following a Hartree-Fock (HF) study of asymmet- 
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ric nuclear matter [15], to be the same as that of the isoscalar part. As a result, quite a 
strong IV mixing was found for 2^ states in ^^'^°0, with the ratio of neutron/proton transi- 
tion matrix elements M^/Mp ~ 4.2 for ^°0. Although that value agrees fairly with previous 
estimates [2, 3, 4] within the limits of experimental errors, a recent measurement of the 
Coulomb excitation of ^°Mg [6] has revealed that the M„/Mp is only around 2.5 for 2^ state 
in ^°0, if one assumes the Mp moment of 2^ state in ^°Mg equal the Mn moment of 2^ state 
in based on the isospin symmetry. A question was raised in Ref. [6] whether such a 
discrepancy is due to the inaccuracy of the (p, p') analysis method of Ref. [12] or the isospin 
impurity in the 2'^ excitation of ^°Mg and ^°0. 

In the mean time, the isovector density dependence of the CDM3Y6 interaction has 
been carefully probed in the CC analysis of the (p, n) reactions exciting the O"*" isobaric 
analog states of ^He [16] and other medium-mass nuclei including ^ospi-^ j-j^'^j^ where the 
isovector coupling was used to explicitly link the isovector part of the nucleon OP to the 
cross section of (p, n) reaction. In particular, a complex isovector density dependence of 
the CDM3Y6 interaction has been constructed based on the microscopic Brueckner-Hatree- 
Fock calculation of nuclear matter [18] by Jeukenne, Lejeune and Mahaux (JLM) before 
being used as folding input. The main conclusion drawn from the results of Refs. [16, 17] is 
that the strength the isovector density dependence of the CDM3Y6 interaction, even after 
it was fine tuned against the JLM results, is somewhat weak compared to the empirical 
isovector strength implied by the (p, n) data. As a result, a renormalization of the (real) 
isovector density dependence of the CDM3Y6 interaction by a factor of 1.2 - 1.3 was found 
[17] necessary to account for the measured (p, n) cross sections. Such an enhancement of the 
isovector density dependence of the CDM3Y6 interaction was also shown [16, 17] necessary 
for a good HF description of the nuclear matter symmetry energy compared to the empirical 
estimates. 

Since a realistic isospin dependence of the effective NN interaction is vital for the deter- 
mination of the IS and IV deformation lengths (or M„ and Mp moments) , a revised folding 
model analysis of the -^^'^'^'^^0(p, j/) data is necessary for a more definitive conclusion on the 
neutron transition strength of 2^ states in 18,20,22q isotopes. After a brief overview of the 
theoretical formalism in Sec. 11, the results of the folding -|- CC analysis of the ^*'^°'^^0(p,p') 
data are presented in Sec. Ill and the main conclusions are given in the Summary. 
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II. GENERAL FORMALISM 



A. Nuclear densities, isoscalar and isovector deformations 

We describe here briefly the method suggested flrst in Ref. [12] to hnk the deformation 
of an excited nucleus and the corresponding transition density based on a collective model 
treatment. As the nuclear deformation is associated with the 'deformed' shape of excited 
nucleus, instead of 'deforming' the optical potential (2), one can build up the proton and neu- 
tron transition densities of a 2'^-pole excitation (A > 2) using the so-called Bohr-Mottelson 
(BM) prescription [19] separately for protons and neutrons 

dp^ (t^ 

Pa(^) = — ' r^p,n. (3) 

Here Pg.s.(r) are the proton and neutron ground state (g.s.) densities and 5^- the corre- 
sponding deformation lengths. Given an appropriate choice of the g.s. proton and neutron 
densities it is natural to represent the IS and IV parts of the total g.s. density as 

p'J:lKr)^plAr)±p',.s.ir). (4) 

One can then generate, using the same BM prescription, the IS and IV parts of the nuclear 
transition density by deforming (4) as 

Pa (0 = -^i) ^ ■ (5) 

The exphcit expressions for the IS and IV deformation lengths are then easily obtained, 

after some integration in parts, as 

N < r^-i >„Sn + Z< r^-i >p 5p 
= A < r^-^ ' 



1 N < r^-i >„ -Z < r^-i >p ' ^ ^ 

Here the radial momenta < r^~^ >x {x — n,p,A) are obtained with the neutron, proton 
and total g.s. densities as 



< r^"^ >x-- 



/oo poo 
plJry^'dr/ plArydr. (8) 

The transition matrix element associated with a given component of nuclear transition 
density (y = n,p, 0, 1) is 

POO 

My= pl{ry+^dr. (9) 



The ratios of the neutron/proton and IS/IV transition matrix elements are given by 



Mp Z < r^-i >p Sp ' 



(10) 



Ml ^ (TV < r^-i >„ -Z < r^-^ >p)di 

Mo {A<r^-^>A)6o ■ ^ ^ 

It is useful to note that there is a one-to-one correspondence between the ratios of transition 
matrix elements in the two representations, and they are related by 

MjMp = (1 + Mi/Mo)/(l - Ml/Mo). (12) 

If one assumes that the excitation is purely isoscalar and the neutron and proton densities 
have the same radial shape (scaled by the ratio N/Z) then Sn = Sp = 5q = Si, 

Mn N Ml N - Z 

TT- — ~^ and — - = — = e. (13) 

Mp Z Mq A ^ ' 

Consequently, any significant deviation of the Mn/Mp ratio from N/Z (or deviation of the 
Mi/Mq ratio from e) would directly indicate an isovector mixing in the nuclear excitation. 

Since the electric transition probabilities B{E2) for 2^ states in i8'20,22q isotopes have 
been measured, we can choose the proton deformation length Sp so that the experimental 
transition rate is reproduced by Be^{E2 f) = e^|Mpp. As a result, the only free parameter to 
be determined from the folding model analysis of the (p, p') data is the neutron deformation 
length Sn- AH the transition matrix elements and other deformation lengths can be directly 
obtained from Sp and Sn using Eqs. (3)-(ll). This feature is the main advantage of our folding 
model approach compared to the standard DWBA or CC analysis using the collective model 
prescription (2). 

We note that the same ^^'^°'^^0(p,p') data have been studied in the folding model using 
the microscopic nuclear transition densities calculated in the quasiparticle random phase 
approximation (QRPA) [20]. In these calculations [3, 5, 21], the QRPA proton transition 
density is scaled to reproduce the experimental B{EX) values, while the strength of the 
neutron transition density is adjusted to the best DWBA or CC fit to the {p.p') data. 
The Mn and Mp transition matrix elements given by the 'scaled' QRPA transition densities 
are then compared with the empirical estimates. Since different effective NN interactions 
were used in the folding calculations of Refs. [3, 5, 21], it is of interest from the reaction 
theory point of view to probe the microscopic QRPA transition densities in our folding model 

6 



analysis using the same effective NN interaction. Therefore, in addition to the BM transition 
densities (3), we have used in the present work also the QRPA transition densities for 2f 
states in Oxygen isotopes given by the continuum QRPA calculation by Khan et al. [20] to 
calculate the inelastic FF. The proton and neutron g.s. densities obtained in the Hartree- 
Fock-Bogoljubov study [22] were used in the folding model calculation of the proton-nucleus 
optical potential. 

B. Folding model with complex CDM3Y6 interaction 

In our version [23] of the folding model, the central nucleon-nucleus potential is evaluated 
in a Hartree-Fock manner as 

^ = ^3\vii\i3 > + < U>Ex|i« >], (14) 

where f d(ex) are the direct and exchange components of the effective A^A^ interaction between 
the incident nucleon i and nucleon j bound in the target A. The antisymmetrization gives 
rise to the exchange term in Eq. (14) which makes the nucleon-nucleus potential nonlocal in 
the coordinate space. To separate the IS and IV contributions, one needs to make explicit 
the spin- and isospin dependence of the (energy- and density dependent) A" A" interaction 

VTi{FX){E, p, s) = Vq^^^\e, p, s) + vf^^^\E, p, s){crcr') 

+v^}^''\e, p, s)(rT') + vf}^''\E, p, s)(cTCT')(Tr'), (15) 

where s is the internucleon distance. The contribution from the spin dependent terms {vio 
and Vii) in Eq. (15) to the central nucleon-nucleus potential (14) is exactly zero for the 
(spin-saturated) Oxygen targets considered in the present work. 

Using a realistic local approximation for the exchange term, the nucleon-nucleus potential 
(14) can be obtained [23] in terms of the isoscalar (t/is) and isovector (t/iv) parts as 

U{E, R) = Uis{E, R) ± Uiy{E, R), (16) 

where the -|- sign pertains to incident neutron and - sign to incident proton. The second 
term in Eq. (16) is the microscopic expression for the Lane potential in Eq. (1) as well as 
its prototype in Eq. (2) for the inelastic scattering FF. Each term in Eq. (16) consists of the 
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corresponding direct and exchange potentials 




(17) 



Ui^{E, R) = j {[p„(r) - pp(r)]<(^, p, s) 
+[p„(H, r) - pp{R, r)]v^^{E, p, s)jo{k{R)s)}d\ 



(18) 



where p{r,r') is the (one-body) density matrix of the target nucleus, with p{r) = p{r,r). 
jo{x) is the zero-order spherical Bessel function and momentum k{R) is determined from 



Here, /i is the nucleon reduced mass, U{R) and Vc{R) arc the nuclear and Coulomb parts 
of the OP, respectively. For a consistent description of the elastic and inelastic nucleon- 
nucleus scattering, one needs to take into account explicitly the multipole decomposition of 
the neutron and proton densities that enter the folding calculation (17)-(18). The details of 
the folding calculation of Uis and Uiy are the same as those given in Ref . [23] , excepting the 
use of a realistic local approximation for the transition density matrix taken from Ref. [24] . 
We note that there exists a more sophisticated version of the single-folding approach, known 
as the g- folding model [25], where the nonlocal exchange potential is treated exactly in 
the Schrodinger equation for the scattered wave, using the explicit wave function for each 
single-particle state \j > taken from the shell model. In this sense, our approach is more 
flexible because one needs to use for the folding input only the total proton and neutron 
densities Prif) and the nuclear densities of any structure model can be used. In particular, 
the use of the 'collective model' prescription (3) for the transition densities has allowed us 
to determine the IS and IV deformations of a nuclear excitation. 

For the effective interaction, we used the density- and isospin dependent CDM3Y6 inter- 
action [14]. While the isoscalar density dependence of the CDM3Y6 interaction has been 
well tested in the folding model analyses of refractive a-nucleus and nucleus-nucleus scat- 
tering (see recent review in Ref. [26]), its isovector density dependence was studied only 
recently in the CC analysis [16, 17] of the charge exchange {p,n) reaction exciting the 0"*" 
isobaric analog states of targets ranging from ^He to ^°^Pb. We recall that the IS density 



k'^{E,R) 



2p 



Cm. 



Re U{R) - Vc{R)]. 



(19) 
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dependence of the CDM3Y6 interaction was introduced [14] as 

C^'^^H^, p, s) = F,s{E, (20) 
Fis{E, p) = Co[l + ao exp(-/3op) - 7op], (21) 

where v^q^^"* (s) are the direct and exchange components of the isoscalar M3Y- Paris inter- 
action [27]. Parameters of Fjs were chosen [14] to reproduce the saturation properties of 
symmetric nuclear matter in the HF calculation. With a linear energy dependence included 
into Co, the IS interaction (20) reproduces very well the empirical energy dependence of the 
IS nucleon OP in nuclear matter [28]. 

For an accurate folding model analysis of the ^^'^°'^^0(p, p') data, it is highly desirable to 
have a complex, density- and isospin dependent NN interaction for the input of the folding 
calculation (16)-(18). Following Ref. [17], we have constructed in the present work, explicitly 
for each energy, an imaginary IS density dependence of the same functional form (21) and 
a complex IV density dependence of the M3Y-Paris interaction 

v^i^''\E, p, .) = Fiv(£;, p)vT''\s), (22) 
Fiv(S, p) = Ci[l + ai exp(-/3ip) - 7ip], (23) 

where the parameters were adjusted to reproduce the JLM density- and isospin dependent 
nucleon OP [18] in the HF calculation of nuclear matter. All radial shapes of v^^^^{s) 
were kept unchanged as derived in terms of three Yukawas from the M3Y-Paris interaction 
[27] (see the explicit expressions for v^Q^^^^is) in Ref. [15]). The isovector part of the folded 
proton-nucleus OP has been used in Ref. [17] as the FF for the (p, n) reaction exciting the 
isobaric analog states, based on the isospin coupling scheme. It turned out [17] that the 
strength of the real isovector interaction (22) is quite weak to account for the observed (p, n) 
data and an enhancement of about 20-30% is needed for a good CC description of the (p, n) 
reaction. Therefore, we have scaled parameter Ci of the real IV density dependence (23) 
by a factor of 1.3 before using for the input of the folding calculation (16)-(18). The final 
parameters of the complex density dependences Fis(tv){E, p) arc presented in Table 1. We 
note that the central proton-nucleus potential (16) is supplemented by the spin-orbital term 
obtained with the folding method of Ref. [23] and the spin-orbital terms of the CDM3Y6 
interaction v]^g^'^{p,s) are assumed to have the same IS and IV density dependences as 
those used for the central terms. All the optical model (OM) and CC calculations have been 



TABLE I: Parameters of the complex IS and IV density dependence of the CDM3Y6 interaction 
defined in Eqs. (21) and (23), respectively. 









■P) 






Im 






E (MeV) 


Co 


ao 


(fm=^) 


70 (fm^) 


Co 


ao 


f3o (fm^) 


70 (fm^) 


24.5 


0.2487 


3.8033 


1.4099 


4.0 


0.1504 


6.0964 


15.503 


4.3931 


43.0 


0.2361 


3.8033 


1.4099 


4.0 


0.0869 


8.9092 


9.9237 


4.2128 


46.6 


0.2336 


3.8033 


1.4099 


4.0 


0.1029 


6.8937 


9.1076 


4.2056 


ReFiv(£;,p) 


Im Fiy{E,p) 


E (MeV) 


Ci 


ai 


(3i (fm3) 


71 (fm-'^) 


Ci 


ai 


Pi (fm3) 


71 (fm'^) 


24.5 


0.2668 


6.3227 


13.725 


-3.8888 


0.2010 


9.6207 


16.053 


-4.3670 


43.0 


0.1490 


9.7964 


10.743 


-4.1147 


0.2315 


6.2846 


13.162 


-4.3612 


46.6 


0.1585 


9.3490 


11.683 


-3.9323 


0.2289 


6.0590 


12.407 


-4.4283 



performed using the CC code ECIS97 written by Raynal [29] . 

III. RESULTS AND DISCUSSIONS 

A. Transition strength of 2^ state in ^^O 

For the Oxygen isotopes under study, the 2^ state (at 1.98 MeV) in ^^0 is the most 
studied one. The core polarization by the two valence neutrons in the 2^ excitation of ^^O 
was shown to be quite strong, with the M„/Mp ratio significantly larger than N/Z. The 
neutron transition strength of 2^ state in ^^O has been measured in several experiments, like 
the (direct) inelastic proton and neutron scattering [8, 30, 31] or inelastic pion scattering 
[32]. The Mn/Mp ratio was often deduced by the (collective model) prescription of Bernstein, 
Brown and Madsen (BBM) [1] which has been checked against the data collected for a 
wide range of single-closed shell nuclei. While the BBM analysis of the low-energy proton 
scattering data seems to favor M^/Mp 1.5 for 2^ state in ^^O [2, 4], the inelastic pion 
scattering data were shown to give a much higher M^/Mp ratio of 2.3 to 2.4 [32]. The 
Mn/Mp ratio deduced from the pion scattering data also agrees fairly with that deduced 
from the measured B[E2) strength of 2^ state in the mirror nucleus ^^Ne, using the isospin 
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symmetry [33] . We note that the DWBA analysis of (p, p') and (n, n') scattering data at 24 
MeV [8] is of particular interest for our study, as it is the only attempt to determine the 
IS and IV deformations for 2'![ state in -"^^O prior to our work [12]. By using prescription 
(2) and assuming 5q to be the average of 5 values given by the {p,p') and {n,n') data, the 
authors of Ref. [8] obtained Sq ~ 1.26 ± 0.06 fm and Si ~ 3.14 ± 1.57 fm which correspond 
to Mn/Mp ft! 1.72 ± 0.70. The coupled-channel effect at the proton energy of 24 MeV was 
shown [8] to affect slightly the deduced deformation parameters. 

To compare our microscopic folding model analysis of ^^0(p,p') data with the collective 
model results of Ref. [8], we have performed a coupled-channel (25f" O^g 3^) analysis 
of ^^0{p,p') data at 24.5 MeV [30] using the (complex) OP and inelastic FF given by the 
folding calculation (16)-(18) for the lowest 2+ and 3^ states in ^^O. By adjusting Mp to the 
experimental B{E2 t) = 45.1 ± 2.0 e^fm^ [34] and B{E3 t) = 1120 ± 11 e^fm^ [35] for the 
first 2+ and 3~ states in ^^O, we obtain 6p — 1.04 ±0.02 and 1.45 ±0.01 fm, respectively, for 
the corresponding proton transition densities (3). To effectively account for the higher-order 
dynamic polarization of the OP by the open nonelastic channels, the (complex) strength of 
the CDM3Y6 interaction is first adjusted to the best CC description of elastic scattering 
data and then is used without any further renormalization to calculate the inelastic FF. 
As a result, the only remaining parameter is the neutron deformation length Sn which is 
determined from the best CC fit to the inelastic scattering data. The renormalization factors 
A^R and Nj of the real and imaginary folded OP were first obtained in the OM analysis of 
the elastic data (see Table II). At 24.5 MeV, the OM fit gives A^r ft 0.91 and A^i ft 0.6. 
These values have changed slightly to A^r ft 0.97 and A^i ft 0.57 when the two-channel 
coupling is taken into account. The folded spin-orbital potential needs a renormalization of 
around 0.5 in both the OM and CC calculations. At higher energy of 43 MeV the best-fit 
Ar factor becomes close to unity while Ai remains around 0.6. One can see from Fig. 1 
that the folded OP gives quite a good description of the measured elastic cross section and 
analyzing power at 24.5 MeV in both the OM and CC schemes. For a comparison, we have 
also performed the OM calculation using the phenomenological OP (parameterized in terms 
of Woods-Saxon potentials) taken from an accurate global systematics [36] by Koning and 
Delaroche (KD). Although KD systematics has been developed for nuclei in the mass range 
24 ^ A ^ 209, our OM analysis shows that it works rather well also for the Oxygen isotopes 
under study. 
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FIG. 1: (Color online) Measured differential cross section and analyzing power of the elastic 
scattering at 24.5 MeV [30] versus the OM results given by the microscopic folded OP and phe- 
nomenological OP taken from the global systematics [36] by Koning and Delaroche (KD). The CC 
results are given by the complex OP and inelastic FF calculated in the folding model. 

The inelastic scattering data at 24.5 MeV [30] for 2^ state in ^^O are compared with 

the results of CC calculation in Fig. 2. Like the earlier folding model study [12], we found in 
the present analysis a significant IV mixing in 2f excitation which leads to Mn/Mp ^ 1.55 
(see Table III). This value is about 25% larger than the ratio implied by the isoscalar limit 
{Mn/Mp = N/Z = 1.25). Compared to the elastic channel, the agreement of the CC results 
with the measured analyzing power of inelastic scattering is rather poor at medium angles. 



12 



TABLE II: Renormalization factors N-^ and Ni of the real and imaginary folded potentials used in 
the OM and CC calculations of elastic and inelastic p+i8'20,22Q scattering. 



OM fit 


CC fit 


Target 


E (MeV) 


Nk 


Ni 


Nls 


(7R (mb) 




Ni 


A^LS 


180 


24.5 


0.91 




0.60 


0.50 


658 


0.97 


0.57 


0.50 


180 


43.0 


1.00 




0.69 


0.50 


545 


1.05 


0.64 


0.50 


20o 


43.0 


1.08 




0.68 


0.50 


573 


1.15 


0.60 


0.50 


220 


46.6 


1.00 




0.72 


0.50 


572 


1.01 


0.70 


0.50 



and that could well be due to a simple treatment the inelastic spin-orbital FF adopted in 
our folding method [23] . However, the inelastic spin-orbital FF does not affect significantly 
the calculated inelastic 2+ scattering cross section which is dominated by contribution from 
the central FF, and the widely accepted procedure is to deduce deformation parameters 
by matching the calculated inelastic scattering cross section to the data. To stress the 
reliability of the folding approach, we have done in parallel the same CC calculation but 
using the collective model form factor (2) determined with the phenomenological OP by 
Koning and Delarochc [36] and the same IS and IV deformation lengths. As expected, the 
form factor (2) was found to strongly underestimate the measured 2+ cross section at large 
angles, in about the same way as established earlier in the folding model studies of inelastic 
heavy- ion scattering [10, 11]. This explains naturally why the IS and IV deformation lengths 
of 2f state in ^^O deduced by Grabmayr et al. [8] in their collective model analysis of the 
same data (5o ~ 1-26 fm and 5i ~ 3.14 fm) are significantly larger than the values deduced 
from our folding model analysis (see Table III). The numerical uncertainties of 6p given in 
Table III are fully determined by those of the measured B[E2) values, while an uncertainty 
of around 5% was assigned to 6n which gives a cross-section shift within the experimental 
errors. The uncertainties of the IS and IV deformation lengths and ratios of transition 
matrix elements were deduced directly from those found for 5p and 5n- 

The neutron deformation length found in the CC analysis of the inelastic p-|-^^0 scattering 
data at 24.5 MeV has been used to calculate the inelastic scattering FF at higher 

energy of 43 MeV. With only strength of the complex CDM3Y6 interaction slightly adjusted 
by the CC fit to elastic scattering data, the folding -|- CC description of the measured (p, p') 
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FIG. 2: (Color online) Measured differential cross section and analyzing power of the inelastic 
scattering at 24.5 MeV [30] versus the CC results given by the microscopic folded FF 
obtained with BM and QRPA transition densities. The collective model result is given by form 
factor (2) obtained with the phenomenological OP by Koning and Delaroche [36] (see more details 
in text). 

data at 43 MeV [3] is quite satisfactory (see Fig. 3) without any further adjustment of 
6n for 2^ state in ^'^O. Given the proton transition strength fixed by the measured B{E2) 
value, we conclude that the neutron transition matrix element M„ should be such that the 
ratio Mn/Mp ^ 1.55 or equivalently Mi/Mq ^ 0.22. The latter is two times larger than 
e = 0.11 and, hence, implies a significant IV mixing in the 2^ excitation. Our present result 
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FIG. 3: (Color online) Elastic and inelastic p+^^O scattering data at 43 MeV [3] in comparison 
with the OM and CC results. Notations for the OP and inelastic FF are the same as used in Figs. 1 
and 2. 

also agrees closely with that given by the BBM analysis of low-energy proton scattering 
data {Mn/Mp ^ 1.5 for 2f state in ^^O) [2, 4]. The microscopic QRPA transition densities 
[20] also give a satisfactory description of the inelastic scattering data under study 

(see Figs. 2 and 3). We note that the continuum QRPA calculation of Ref. [20] strongly 
underestimates the E2 strength and gives B{E2 |) = 14 e^fm^ for 2f state in ^^O compared 
to the experimental value of 45 e^fm^. On the contrary, the predicted neutron transition 
strength is much too high [Mn/Mp ^ 2.88) compared to that found in the present work and 
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TABLE III: Deformation lengths and ratios of the transition matrix elements for 2^" states in 



18,20,22 o deduced from the present folding + CC analysis of inelastic proton scattering [see defini- 
tions in Eqs. (3)-(ll)]. 



Nucleus 


N/Z 


e 


5p (fm) 


5n (fm) 


Mn/Mp 


^0 (fm) 


6i (fm) 


Ml /Mo 


180 


1.25 


0.11 


1.04 ±0.02 


1.23 ±0.06 


1.55 ±0.08 


1.15 ±0.04 


1.87 ±0.27 


0.22 ±0.03 


20o 


1.50 


0.20 


0.82 ± 0.03 


1.95 ±0.10 


3.24 ±0.20 


1.52 ± 0.06 


3.75 ± 0.26 


0.59 ± 0.05 


220 


1.75 


0.27 


0.70 ± 0.12 


0.91 ± 0.05 


1.81 ±0.33 


0.85 ± 0.05 


1.14 ±0.16 


0.43 ± 0.07 



other studies. Such a strong neutron transition strength seems to compensate for the weak 
proton transition strength predicted by the QRPA, and the inelastic FF folded with the 
QRPA transition densities gives a reasonable description of the inelastic scattering data. It 
is interesting that about the same good description of the inelastic scattering data is given 
by the renormalized QRPA densities, with p^(r) scaled to reproduce the measured B{E2) 
value and P2{r) scaled to give the same M„ as that obtained above for the BM model (3) of 
transition densities (see Fig. 2 and lower panel of Fig. 3). 

As mentioned above, an alternative method to determine the neutron transition matrix 
element M„ has been suggested some 30 years ago by Bernstein et al. [33] based on the 
isospin symmetry. Namely, M„ can be obtained from Mp measured for the same 2^ excitation 
in the mirror nucleus with an electromagnetic probe if one assumes the charge independence 
of the 2+ excitation in members of a T-isospin multiplet. In particular, one has for the 
isobars with opposite signs of the isospin projection 

Mp{-T,)^Mn{T,). (24) 

Using the electric transition rates B{E2) for ^^O and ^^Ne taken from the latest compila- 
tion of the experimental data [34] , we easily deduce the ratio of transition matrix elements 
Mn/Mp fti 2.44 ± 0.18 for 2f states in ^^O using Eq. (24). This value is significantly larger 
than that obtained in the present folding + CC analysis and collective model analyses re- 
ported in Refs. [2, 4, 8]. To illustrate such a difference in terms of the {p,p') cross section, 
we have done the same folding + CC calculation but using the enhanced neutron transition 
density (3) which gives Mn/Mp 2.44. Then, the calculated inelastic cross sections strongly 
overestimate the measured (p, p') data at both energies under study (see Fig. 4) . Assuming 
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FIG. 4: (Color online) Inelastic scattering data at 24.5 [30] and 43 MeV [3] for 2f state 

in ^^O in comparison with the CC results. The neutron deformation length 6n of the transition 
density (3) was adjusted to give two different ratios Mn/Mp = 1.55 and 2.44. 

the realistic value Mn/Mp ^ 1.55, we might interpret the difference shown in Fig. 4 as an 
indication to the isospin impurity in the 2^ excitations of the A = 18, T = 1 isobaric mul- 
tiplet. It is complementary to note that a similar isospin- impurity effect has been found by 
Khan et al. [37] for the 2^ excitations of the A = 30, T = 1 isobaric multiplet. Neverthe- 
less, if one takes into account the data of inelastic vr"'' and vr^ scattering from ^^O [32] which 
give Mn/Mp ^ 2.3 — 2.4 in a distorted-wave impulse approximation analysis using different 
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types of the transition density for 2'1 state, then the measured B{E2) value of 2^ state in 
the mirror ^^Ne nucleus seems to support a good isospin symmetry in this case. Given the 
accurately measured electric transition rates B{E2) of 2^ states in ^^O and -"^^Ne, a future 
(high-precision) experiment to re-determine the neutron transition strengths of 2^ states 
in these two mirror nuclei should provide vital data for the determination of the isospin 
impurity using relation (24). 

B. Transition strength of 2'1 state in 

In difference from the stable ^^O target, inelastic proton scattering to 2f state in the 
unstable isotope has been measured only recently, in the inverse kinematics, at 30 MeV 
[2] and 43 MeV [3]. Using the measured transition rate B{E2 ]) ^ 28.1 ± 2.0 e^fm^ [34] 
for 2]^ state in ^°0, we have deduced the proton deformation length 5p ~ 0.82 ± 0.03 fm 
for the proton transition density (3). The proton deformation length of 3]^ state in was 
taken from the empirical estimate of Ref . [12] , and the neutron deformation length was then 
adjusted to the best CC fit to the inelastic ^^0(p,p') data for 3f excitation [3]. Similar to 
our earlier folding -|- DWBA analysis [12] of these data, the best-fit M„/Mp ratio for the 3^ 
excitation turned out to be quite close to N/Z = 1.5 which indicates a dominant IS character 
of this state. The isospin character of 2]*" state is very much different from that found for 
3]^ state. Using the deformation parameters extracted from the CC analysis of the 30 MeV 
data with the collective model FF, the ratio of neutron/proton transition matrix elements 
M„/Mp 2.9±0.4 has been deduced by Jewell et al. [2] for 2t state in ^OQ. The JLM folding 
model analysis of 43 MeV data by Khan et al. [3] using the QRPA transition densities has 
found a stronger IV mixing in this state, with Af„/Mp ^ 3.25 ± 0.80. Our earlier folding 
model analysis [12], using the compact method (3)-(ll) and original CDM3Y6 interaction 
with the IV density dependence assumed to be the same as the IS one, has given a larger 
ratio of M^/Mp ~ 4.2 for 2'^ states in ^°0. In the present work we concentrate on the 43 
MeV data which contain more data points and cover a wider angular range. The results 
of our folding -|- CC analysis are compared with the elastic and inelastic p+'^^O scattering 
data at 43 MeV in Fig. 5. With 6p fixed above by the measured B{E2) value, the best-fit 
neutron deformation length is 5„ ~ 1.95 fm which results on the ratios M„/Mp ^ 3.24 and 
Ml /Mo w 0.59. These values are well exceeding the IS limit of M^/Mp fa N/Z =1.5 and 
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FIG. 5: (Color online) Elastic and inelastic p-\-'^^0 scattering data at 43 MeV [3] in comparison 
with the OM and CC results. Notations for the OP and inelastic FF are the same as used in Figs. 1 
and 2. 

Mi/Mq ~ e = 0.2. The deduced IV deformation length (see Table III) is about 2.5 times 
the IS deformation length and confirms, therefore, a strong core polarization by the valence 
neutrons in 2^ excitation of the open-shell ^^O nucleus. In difference from the inelastic 
scattering data at 24.5 MeV shown in Fig. 2, the 43 MeV data for i8,20q 'targets' 
(measured at angles 6c.m. < 50° only) are reasonably reproduced by both the folded and 
collective model form factors based on the same 6o and 6i. A substantial difference between 
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the CC results given by these two choices of inelastic FF was found at larger scattering angles 
(50° < ©c.m. < 180°) where no data point was taken. In this sense, a future experiment 
aiming to measure {p,p') cross section over a wider angular range could provide a better 
test ground for the inelastic FF and neutron transition strength. 

The continuum QRPA description of 2f state in ^'^O is better than that for ^^O, with the 
predicted B{E2 f) ~ 22 e^fm^ (compared to the adopted experimental value of 28 e^fm^) 
and Mn/Mp 3.36 [20]. A good CC description of the (p,p') cross section was obtained 
after a slight renormalization of the QRPA transition densities to reproduce the experimental 
B{E2) value and best-fit Mn/Mp ratio (see Fig. 5). We stress that the use of a more reahstic 
version of the (complex) density- and isospin dependent CDM3Y6 interaction in the present 
work has pinned down the best-fit ratio of transition matrix elements for 2^ state in to 
Mn/Mp 3.24 which is very close to that deduced from the JLM folding model analysis [3] 
of the same data. Together with the results of the BBM analysis reported in Refs. [2, 4], 
our results confirm again a strong IV mixing in 2^ excitation of ^'^O. 

Like the ^^O case, there is an alternative method to determine the Mn/Mp ratio from the 
electric B(E2) transition rates measured for 2]*" states in ^^0 and its mirror partner ^°Mg, 
using Eq. (24) given by the isospin symmetry. Towards this goal, a measurement of the 
Coulomb excitation of the unstable ^°Mg nucleus has been performed by RIKEN group [6] 
using a radioactive ^°Mg beam incident on the lead target. The extracted transition rate 
B{E2 t) 177 ± 32 e^fm^ for 2^ state in ^'^Mg seems to agree well with the prediction of 
realistic cluster model for this nucleus. Assuming the proton transition matrix element Mp 
for 2f state in ^°Mg equal the neutron transition matrix element M„ for 2f state in ^°0, 
one obtains easily Mn/Mp ~ 2.51 ± 0.25 for the latter. This value is about 30% smaller 
than the best-fit Mn/Mp ratio obtained in Ref. [3] and present work. Given an accurate 
treatment of the folding model ingredients and similar effect discussed above for 2'^ state 
in ^^O, such a difference in Mn/Mp ratios deduced by the two methods might well indicate 
the isospin impurity in the 2^" excitation of the A = 20, T — 2 isobaric multiplet. To 
illustrate this effect in the calculated inelastic cross sections, we have plotted in Fig. 6 
the CC result obtained with two choices of the neutron transition density (3) which were 
scaled to give Mn/Mp ^ 2.51 and 3.24. One can see that the CC results associated with 
Mn/Mp 2.51 substantially underestimate the measured {p,p') data. In terms of the 
total (p,p') cross section, the difference caused by the 'isospin impurity' shown in Fig. 6 
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FIG. 6: (Color online) Inelastic p+^^O scattering data at 43 MeV [3] for 2f state in ^^O in 
comparison with the CC results. The neutron deformation length 6n of the transition density (3) 
was adjusted to give two different ratios Mn/Mp = 2.51 and 3.24. 

is around 40%. This difference reduces to around 30% when one adopts the upper limit 
of the measured transition rate, B{E2 |) ^ 210 e^fm"^, for 2^ state in ^''Mg [6]. We must 
note, however, that the last 2 data points in Fig. 6 seem to agree better with the CC results 
associated with Mn/Mp ^ 2.51. In the same logic as discussed above for ^^O and -"^^Ne, a 
future inverse-kinematics measurement of ^°0(p, p') and ^°Mg(p, p') reactions to determine 
neutron transition strength of 2^ states in these two mirror unstable nuclei could provide 
important data for the check of isospin purity in the 2f excitation of the A = 20, T = 2 
isobars using relation (24). 
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C. Transition strength of 2 J" state in ^^O 

If we consider as consisting of the ^^O core and four valence neutrons, then the 
large IV deformation length Si extracted above for 2'^ state in indicates a strong core 
polarization by the valence neutrons in the 2^ excitation. In such a 'core + valence neutrons' 
picture, it is natural to expect that 2f state in ^^O should be more collective and have a 
larger IV deformation length due to the contribution of two more valence neutrons. However, 
the inelastic ^^0{p,p') scattering data at 46.6 MeV measured recently at GANIL [5] show 
clearly the opposite effect, with the inelastic cross section about 3 to 4 times smaller than 
^^0(p,p') cross section measured at 43 MeV [3] for 2^ state in ^OQ. The folding + DWBA 
analysis [5] of these data using the method (3)-(ll) and original CDM3Y6 interaction (with 
the IV density dependence assumed to be the same as the IS one) has pointed to a much 
weaker neutron transition strength of 2f in ^^O. Given a significantly higher excitation 
energy of this state (1.5 MeV higher than that of 2f state in ^°0), the ^^0(p,p') data at 46.6 
MeV were considered [5] as an important evidence for the neutron shell closure occurring 
at = 14 or 16. Given the measured transition rate B{E2 t) ~ 21 ± 8 e^fm^ [34] for 2^ 
state in ^^O, we have obtained the proton deformation length Sp 0.70 ± 0.12 fm for the 
proton transition density (3). Like in cases of ^^'^°0, the neutron deformation length for 
2^ state in was adjusted to the best description of the 46.6 MeV data by the folding 
+ CC analysis. Since there are no data measured for {p,p') scattering at 46.6 MeV to 3]~ 
state in ^^O, we have considered in our CC scheme only the coupling between the elastic and 
2+ inelastic scattering channels (0+g <-> 2^). The folding + CC results are compared with 
the elastic and inelastic scattering data in Fig. 7. With Sp fixed by the measured 

B{E2) value, the best-fit neutron deformation length was found 5n ~ 0.9 fm which lead 
to the ratios Mn/Mp ^ 1.81 and Mi/Mq ^ 0.43. These values are rather close to those 
impfied by the IS hmit of Mn/Mp ^ N/Z = 1.75 and Mi/Mq ^ e = 0.27. The deduced IV 
deformation length Si is around 30% larger than the IS deformation length Sq (see Table III) 
and this is much smaller than the difference between Si and Sq found above for 2'1 state in 
^°0. Thus, our results show a much weaker polarization effect by the valence neutrons in 
the 2f excitation of ^^O nucleus. This subtle effect could not be accurately described by 
the continuum QRPA. While the QRPA calculation gives B(E2 |) f« 22 e^fm"^ (in a perfect 
agreement with the measured value [34]), the predicted Mn/Mp 3.53 [20] is nearly two 
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FIG. 7: (Color online) Elastic and inelastic p+'^^O scattering data at 46.6 MeV [5] in comparison 
with the OM and CC results. Notations for the OP and inelastic FF are the same as used in Figs. 1 
and 2. 

times the empirical data. That is the reason why the CC results given by the original QRPA 
transition densities for 2^ state in ^^O strongly overestimate the data as shown in Fig. 7. 
Like the case, the inelastic data are reasonably reproduced by both the folded 

and collective model FF based on the same deformation IS and IV lengths. A significant 
difference given by the two choices of inelastic FF can be seen at larger scattering angles 
where no data point was measured. 
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FIG. 8: (Color online) Inelastic p+^'^'^^O scattering data versus the CC results given by inelastic 
folded FF with the isospin dependence of the CDM3Y6 interaction included (IS+IV) or neglected 
(IS only). 

Given a smaller IV deformation resulted from a weaker core polarization by the valence 
neutrons in 2^ excitation of ^^O compared to the case, it is quite illustrative to show the 
explicit IV contribution in the calculated inelastic cross section. We have plotted in Fig. 8 the 
calculated inelastic p+^*''^^0 cross sections given by inelastic folded FF with the contribution 
from the isovector part (18) included or neglected. One can see that the IV contribution is 
very strong and amounts up 40% of the total inelastic cross section in the case, and it 
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becomes weaker (around 22%) in the ^^O case. As can be seen from Eq. (18), the IV part 
of the folded FF is entirely determined by the difference between the neutron and proton 
transition densities and its strength is, therefore, directly proportional to the contribution by 
the valence neutrons. Consequently, it is vital to treat the IV dependence of the (effective) 
NN interaction properly in a folding model analysis of proton scattering on a neutron-rich 
target before the neutron transition strength can be accurately deduced. Although one 
could still obtain a good description of the (p,p') data with the IS form factor only by 
scahng up the strength of transition densities, hke in the folding model analysis [21] of these 
same data using the same continuum QRPA transition densities and isospin independent 
DDM3Y interaction, it is uncertain to compare the best-fit Mn/Mp ratios obtained in such 
an analysis with those deduced by our consistent folding method. We note further that 
the more advanced g-folding model [25] was also used recently to study the same inelastic 
proton scattering data on Oxygens [38]. While the inelastic cross sections given by this 
g-folding model agree fairly with data in the angular shape, the authors need to scale the 
calculated (p,p') cross sections by a factor of 2, 5 and 1.6 to fit the data for 2f excitation 
in ^^O, and ^^O, respectively, and this is likely due to a truncated single particle basis 
based on Ohu = shell model wave functions only. The fact that the largest scaling was 
needed for the case in the g-folding model study also indicates that the higher order 
configuration mixing caused by the valence neutrons is strongest in and this result is in 
a sound agreement with our finding. Finally, it is interesting to note that the contribution 
by the IV form factor to the (p,p') cross section for 2f excitation of is very close to 
the cross section shift presumably caused the 'isospin impurity' in the 2f excitation of the 
A = 20, T = 2 isobars (compare Fig. 6 and upper panel of Fig. 8). 

With a direct connection between the IV deformation and dynamic contribution by the 
valence neutrons to the nuclear excitation, it is natural to link the IV deformation with 
possible changes of the neutron shell structure. The best-fit IS and IV deformation lengths 
of 2f states in Oxygen isotopes and those derived from the results of continuum QRPA 
calculation [20], using Eqs. (9) and (10), are plotted versus the neutron number N in Fig. 9. 
For the double-closed shell ^^O nucleus, we have adopted the IS limit with 6q = 6i = 
1.038 ± 0.048 fm as deduced from the measured B[E2) value and used in a recent folding 
model study of inelastic ^^0-1-^^0 scattering [39]. An enhanced IV deformation (with 5i > 5q) 
resulted from the core polarization by the valence neutrons can be seen for the open-shell 
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FIG. 9: (Color online) Observed excitation energies (upper panel), the isoscalar (Jq) and isovector 
{6i) deformation lengths (lower panel) of 2'^ states in Oxygen isotopes deduced from the folding 
model analyses of this work and Ref. [39] (for double-closed shell ^^O) and from the continuum 
QRPA results [20]. 

i8,20q n^cigi^ with maximum of 6i observed for 2^ state in or at = 12. Such a 
maximum of the IV deformation also corresponds to the largest M„/Mp ratio found for 
2f state in ^°0. With N approaching 14, the extracted 6i value is drastically reduced and 
becomes rather close to Sq which indicates a much weaker IV mixing in 2^ excitation of ^^O. 
A similar trend has also been predicted by the continuum QRPA calculation [20], although 
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the predicted difference between 5i and 5o still remains significant at N — 14. This difference 
was predicted to be substantially smaller at N — 16 (see open squares and circles in Fig. 9), 
and it is natural to suggest from the A^-dependence shown in Fig. 9 that Si is reaching 
its second minimum at = 16. Based on a similar N dependence obtained for the IV 
deformation lengths of 2^ states in Sulfur isotopes [13] where a clear minimum of 6i was 
found at the neutron magic number = 20, the deduced N dependence of Si for 2f states 
in Oxygen isotopes seems to suggest that the neutron shell closure occurs again at A?" = 16 
in Oxygen isotopes. Such a shell closure scenario is well illustrated by the N dependence 
of the excitation energies of 2i states in Oxygen isotopes (upper panel of Fig. 9). One 
can see that the energy of 2^ state goes through its minimum at A^ = 12, where the core 
polarizing contribution by the valence neutrons is strongest. As A^ moves to A^ = 16 this 
contribution becomes much weaker and the excitation energy of 2f state becomes larger 
because of the enhanced energy gap between the Is and Od neutron subshells. Although the 
excitation energy of 2f state in ^^O was predicted by different structure calculations [20, 40] 
to be around 4 MeV, the experimental observation has been quite difficult due to the weak 
excitation of this state. In particular, no (p,p') data could be measured so far for 2f state in 
^''O. An important evidence has been found recently in the experiment on neutron decay of 
unstable Oxygen isotopes by the Michigan State University Group [41, 42], where 2f state 
of ^^O was identified as a very weak resonance at the excitation energy of about 4.7 MeV 
which undergoes direct neutron decay to the ground state of ^^O. Assuming the weak E2 
transition strength predicted by the continuum QRPA [20] for 2^ state in ^^O, our folding 
model approach predicts that the {p,p') cross section for this state is at least factor of 2 
smaller than that measured for 2^ state in ^^O [5]. The (unbound) excited 5/2+ state of ^^O 
has also been observed in the same neutron decay measurement [41] at an excitation energy 
of around 2.8 MeV which fits well into the gap of about 4 MeV between the lsi/2 and Ods/2 
subshells predicted, e.g., by the Hartree-Fock-Bogoljubov calculation [40]. Based on this 
discussion as well as the systematics on the /3-decay Q values and single neutron separation 
energies made by Kanungo et al. [43] for a wide range of neutron rich even-even isotopes, we 
can draw a definitive conclusion on the neutron shell closure at A^ = 16 in unstable Oxygen 
isotopes. To this end, more experiments for ^^O, especially, the {p,p') measurement in the 
inverse kinematics would be of further interest. 
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IV. SUMMARY 



A coupled channel analysis of the ^^'^°'^^0(p,p') scattering data has been performed, using 
the OP and inelastic FF calculated microscopically in a compact folding model approach, to 
extract the neutron transition matrix elements M„ as well as the isoscalar (Sq) and isovector 
(Si) deformation lengths of 2f states in the Oxygen isotopes, with the proton transition 
matrix elements Mp fixed by the measured electric transition rates B[E2). 

The newly determined ratios Mn/Mp for 2]*" states in ^^'^OQ have been compared to 
those deduced from the isospin symmetry, using the experimental B{E2) transition rates 
of 2'1 states in the proton rich ^^Ne and ^°Mg isotopes. Given the experimental B{E2) 
values available for 2'^ states in these four nuclei, a future high- precision (p, p') measurement 
to accurately determine the neutron transition strengths of 2f states in the mirror pairs 
i8Q^i8j^g and ^°0,^''Mg (in the inverse kinematics), should provide vital information on the 
isospin impurity in the 2^ excitation of ^4 = 18, T = 1 and A = 20,T = 2 isobaric multiplets, 
respectively. 

The enhancement of the IV deformation has been confirmed again for the open-shell 
i8,20q nuclei -vvhich show a strong core polarization by the valence neutrons. Along the 
isotope chain, the behavior of the dynamic IV deformation of 2^*" state is closely correlated 
with the evolution of the valence neutron shell, and 6i has been found to reach its maximum 
at = 12 which corresponds to the largest Mn/Mp ratio found for 2f state in ^°0. A fast 
decrease of the IV deformation towards A^ = 16 should be connected with the neutron shell 
closure occurring at this new magic number of neutrons. 
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